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Abstract 


A  Slaulatlon  Gass  to  train  data  processors  In  optlaal  data 
recovery  techniques  la  Included  In  this  report.  Also,  «e  studied 
the  trlple-4ultlplexlns  respons*  to  a  point  source  for  a  future 
scanning  spectroaetcr.  Three  basic  notions  are  required  to  nodu- 
late  the  wavelength  and  angular  coordinates  of  the  source.  When 
pesfomlng  prellnlnary  cxperlnents  with  the  Multiplex  Scanner, 
we  discovered  that  the  underlying  Talbot  effect  la  useful  also  for 
neasurlng  phase  objects.  The  theory  of  this  new  "Talbot  Interfere 
oneter"  Is  explored  and  experlaental  results  of  a  candle  flans 
are  reported.  For  axially  aynnetrlc  objacta,  thla  setup  Is  nodl- 
fled  by  replacing  the  straight  line  gratings  with  circular  gratings. 
In  this  way  the  radial  gradient  of  the  deformed  wavefront  la 
observed.  Finally  a  naans  of  simply,  yat  accurately  colllnatlng 
a  1 Ight  source  using  the  Talbot  Interferometer  Is  described. 
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Introduction 


The  goal  is  to  measure  infrared  radiation  from  a  balloon  at 

an  elevation  of  30  km.  Maximum  coverage  of  the  sky  by  the  scanning  spectro¬ 
meter  during  the  flight  of  the  balloon  is  desired.  To  this  end  our  previous 
Annual  Report  of  30  June  1970  devoted  its  efforts  in  describing:  (1)  the 
optimum  scanning  mode  for  the  existing  LMSC  instrument*  (2)  a  scheme  for  pro¬ 
cessing  the  data  gathered  by  this  instrument  and  (3)  a  future  instrument  based 
on  the  Talbot  effect.  Our  efforts  since  then  have  been  to  support  data  pro¬ 
cessing  for  the  existing  equipment  (see  S2)  and  to  continue  the  exploration  of 
future  methods  of  Instrumentation  (§1). 

During  our  investigation  of  the  scanning  spectrometer  we  found  that  the 
collimating  lens  had  to  be  good.  Furthermore,  it  was  discovered  that  the  Talbot 
interferometer  cculd  be  set  up  to  test  the  lens  and  other  optical  components 
(S3  and  Appendix  A).  The  straight  line  gratings  of  our  setup  laterally  shears 
the  wavefront  that  passes  through  the  test  object,  giving  multiple- interferences. 
Our  theory  showed  that  if  we  allow  only  the  zeroth  diffracted  order  to  pass  to 
the  plane  of  observation,  the  second  derivative  of  the  object  transmittance 
would  be  observed.  On  the  other  hand  if  either  of  the  first  diffracted  orders 
were  allowed  to  pass,  then  we  would  observe  the  first  derivative.  Good  colll- 
mation  of  the  light  source  was  a  critical  factor  in  this  experiment.  By 
separating  the  two  gratings  to  a  maximum  distance  over  which  interference  could 
still  be  observed,  a  very  sensitive  test  of  the  degree  of  colllmatlon  was 
obtained  (§4  and  Appendix  B).  Objects  with  axial  syasMtry  are  best  tested  by 
interferometers  that  radially  shear  the  deformed  wavefront.  For  this  reason 
we  modified  our  Talbot  Interferometer  by  replacing  the  straight  line  gratings 
with  circular  gratings.  Our  theory  and  experiments  are  not  quite  complete. 

We  report  here  on  the  theory  of  the  Talbot  effect  in  cylindrical  coordinated 
for  axially  symmetric  objects  (§5). 


§1  The  Scanning  Spectrometer 


In  the  previous  Annual  Report  of  30  June  1970,  a  spectroiseter  based  on  the 
Talbot  effect  was  described.  A  schematic  representation  of  this  Instniment  is 
shown  in  Figure  1.  Sources  to  be  observed  are  assumed  to  be  sufficiently  far 
away  so  that  waves  Impinging  upon  the  first  grating  are  essentially  plane.  We 
use  the  colllsmtlng  lens  and  point  source  to  simulate  a  sky  scene.  The  plane 
wave  is  diffracted  by  the  first  grating  and  |;he  Talbot  images  so  formed  are 
detected  by  the  second  grating.  The  photodiode  receiver  collects  the  light 
which  passes  through  both  gratings.  There  are  three  variables  of  Interest  in 
the  spectrum  of  a  source,  namely  the  wavelength  and  the  elevation  and  azimuth 
angles  that  locate  the  direction  from  which  the  radiation  is  emitted.  Since 
we  have  essentially  only  one  signal  parameter  (i.e.  the  intensity  of  light 
source  on  the  photodiode  current),  a  triple-multiplexing  scheme  is  required  in 
order  to  place  the  three  parameter  function  onto  the  one-parameter  signal.  The 
following  motions  of  the  spectrometer  will  accomplish  this  task.  The  movement 
of  the  second  grating  in  the  x-directlon  is  used  for  obtaining  cosinusoidal 
coding  of  the  frequency  of  the  spectrum.  Rotation  of  the  Instnaent  will  encode 
the  angular  position  of  the  sources  (a,  S)  and  z-motlon  of  the  second  grating 
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will  Increase  the  delay  time  for  good  apectral  resolution  and  la  necessary  to 
recover  the  spectrum. 

Our  studies  of  the  triple-multiplexing  response  to  a  point  source  for 
this  Instrument  show  that  we  have  three  sialn  carrier  frequencies  produced  as  a 
consequence  of  the  three  Independent  motions.  They  are:  "  ^^x/d;  <^long 

■  irv^A/d^;  and  Uang  *  'P-  Here  v^  and  V2  are  the  velocities  In  the  x-  and  z- 
dlrectlons,  d  Is  tne  grating  period,  X  is  the  wavelength  of  the  radiation,  Is 
the  time  derivative  of  the  angular  position  q>  of  the  point  source  as  shown  In 
Fig.  1,  and  oi^Qg  are  the  carrier  frequencies  associated  with  the 

y-  and  z-dlrection  and  the  rotational  motions  respectively.  The  relative 
uotlons  will  have  to  be  regulated  In  a  manner  so  that  the  three  main  carrier 
frequencies  can  be  separated.  This  In  turn  will  be  determined  by  the  desired 
resolutions  of  the  three  spectrum  parameters  (X,  a,  B).  The  analysis  Is  Incom¬ 
plete  but  a  detailed  report  will  be  made  later. 


§2  A  Simulation  Game 

A  "Signal  Processing  Game"  was  devised  for  the  purpose  of  training  the 
person  In  charge  of  signal  processing  whom  we  call  the  "Investigator".  The 
game  consists  of  digital  data  (64  samples)  which  describe  a  Signal  plus  noise. 
The  Investigator  Is  supposed  to  devise  a  strategy  for  Impvovlng  the  slgnal-to- 
nolse  ratio.  In  problem  #1  the  Investigator  knows  very  little  about  the 
statistical  behavior  of  signal  and  noise.  Hence  only  a  mild  lsq>rovament  Is 
possible.  In  the  subsequent  problems  the  Investigator  Is  supplied  with  more 
and  more  statistical  Information,  which  allow  for  better  and  more  sophisticated 
data  Improvement  schemes. 

The  rules  request  that  the  game  keeper  hand  the  "Introduction"  and 
"Recorded  Data"  sections  to  the  Investigator.  Next  the  game  keeper  gives  the 
Investigator  problem  #1.  It  Is  Important  that  only  one  problem  be  given  the 
Investigator  at  a  time,  because  the  forsnilatlon  of  problem  #2  reveals  part  of 
the  answer  to  #1.  Also  the  last  section,  "About  the  Signal  Synthesis",  should 
be  well-hidden  from  the  Investigator  since  It  explains  everything  which  the 
"Investigator"  shall  find  for  himself,  aatong  others  the  correct  noise-free 
signal. 
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S2.1  Introduction  to  a  Signal  Proccaslng  Gaae 


The  ale  of  this  gaae  la  to  develop  akllls  lil  signal  processing.  At  the 
saae  tlae  errors  in  the  computer  program  for  this  Job  will  become  apparent. 

The  definitions  are  essentially  those  of  the  1970  final  report.  Notice: 
there  is  an  error  on  pg.  38  of  that  report.  In  the  denominator  of  equ.  (13.18) 
|ng|^  should  read  Ihg/EP. 

The  input  data  for  this  game  are  the  (recorded  data  U|^(t).  It  is  assumed 
that  the  non-uniform  scanning  velocity  of  the  spinning  gondola  has  been 
corrected  for  already.  The  time  coordinate  is  given  in  discreet  numbers  t  ■ 

0,  1,  2... 63.  We  may  consider  Ug(t}  as  being  about  one  quarter  of  a  single 
horizontal  scan  (8  -  constant) . 

The  rules  of  the  games  are  as  follows.  The  "investigator"  gets  the 
sheet  "Recorded  Data  Ug(t)"  and  the  sheet  "Problem  #1".  After  solving  this 
problem  he  will  give  the  solution  to  the  "monitor"  and  to  the  "game  constructor". 
Nov  he  Bay  start  on  problem  #2,  and  so  on.  But  it  is  vnportant  that  the  inves¬ 
tigator  does  not  get  the  next  problem  sheet  before  he  has  finished  the  previous 
problem.  The  reason  is  that  the  formulation  of  the  later  problems  contains 
parts  of  the  answers  to  the  earlier  problems.  This  has  to  do  with  the  basic 
structure  of  this  simulation  game:  for  performing  any  meaningful  signal  pro¬ 
cessing  operation  one  must  have  some  knowledge  about  the  original  signal  and/or 
the  noise.  For  example  in  problem  #1  the  investigator  is  told  that  the  noise 
is  additive  and  non-negative.  In  the  later  problems  the  investigator  will  be 
supplied  with  even  more  a  pvnovi  Information.  Naturally,  this  should  enable 
him  to  extract  the  signals  better  and  better.  But  the  methods  for  doing  this 
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Increase  in  coaplexlty. 


On  the  very  last  pages,  following  the  probleas.  the  design  of  the 
"recorded  data"  is  explained,  and  the  true  original  signal  is  unveiled. 
Obviously  those  pages  should  not  be  given  to  the  investigator  before  he  has 
solved  all  the  probleas. 


§2.2  Recorded  Data  U»(t) 

t  is  the  discrete  tlae  variable,  running  froa  t«0  to  t*63. 


t 

u.(t) 

t 

Ug(t) 

t 

u.(t) 

t 

0 

55 

16 

110 

32 

39 

48 

383 

1 

25 

17 

184 

33 

56 

49 

10 

2 

85 

18 

29 

34 

05 

50 

69 

3 

61 

19 

51 

35 

15 

51 

58 

4 

20 

20 

42 

36 

95 

52 

52 

5 

95 

21 

78 

37 

09 

53 

66 

6 

07 

22 

09 

38 

81 

54 

79 

7 

00 

23 

15 

39 

21 

55 

134 

8 

62 

24 

13 

40 

81 

56 

94 

9 

79 

25 

50 

41 

399 

57 

102 

10 

148 

26 

99 

42 

372 

58 

108 

11 

105 

27 

54 

43 

348 

59 

94 

12 

125 

28 

99 

44 

303 

60 

56 

13 

125 

29 

35 

45 

383 

61 

67 

14 

173 

30 

98 

46 

317 

62 

51 

15 

181 

31 

02 

47 

317 

63 

63 

S2.3  Problems 


Problem  #1 

Given  are  the  recorded  data  UK(t)  with  t  ■  0,  1...63.  Wanted  are  the 
original  data  Uo(t),  which  represent  the  “one-dlmenslonal  equivalent  object 
radiation"  Sg(a) .  We  assume  that  the  known  Influences  of  the  telescope 

R(x',y')]  and  of  the  electrical  system  [G(t)]  have  been  compensated 
already  or  are  negligible.  But  the  recorded  signal  Uj((t)  Is  corrupted 
by  additive  noise  N(t): 

UR(t)  -  Uo(t)+N(t). 

The  only  features  known  about  the  original  signal  Uo(t)  and  about  the  noise 
N(t)  are  that  they  are  non-negative: 

Uo(t)  >  0;  N(t)  i  0. 

Furthermore  the  noise  N(t)  Is  stationary,  which  means  that  the  noise  properties 
are  not  "drifting".  In  other  words,  short-term  average  features  of  the  noise 
remain  the  same  from  the  beginning  to  the  end  of  the  observation. 

Try  to  utilize  the  given  a  priori  Information  for  computing  a  new  signal 
Ui(t)  from  UR(t),  which  somehow  Is  better  then  U^Ct)  es  an  approximated 
representation  of  Uo(t).  Plot  U^(t)  as  a  continuous  curve,  and  also  U^Ct)  for 
comparison. 
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Problea  #2 


Given  are  the  facts: 

Uo(t)  >  0;  N(t)  >0;  N  -  50. 

By  N  we  aean  the  linear  average  of  the  noise.  This  N  can  be  visualized  as  the 
dark  current  of  the  photoreceiver  as  measured  with  an  Instrument  which  rejects 
high  frequencies. 

Baaed  on  these  facts  try  to  compute  a  better  signal  U^Ct)  from  U^Ct). 
Plot  both  1)2  (t)  and  Ug(t) . 


(to  be  cut  by  game  monitor) 


Problem  #3 


Given  are  the  same  facts  as  In  the  previous  problem.  In  addition  It  Is 
known  that  the  noise  Is  approximately  "white”. 


N(t)  -  N  +  h(t);  n(v) 


■ 


n(t)  e 


-2irlvt 


dv; 


V  -  m/64;  m  -  -32,  -31, . . .-1,0,+1, . . .+30,  +31; 

|fl(v)P  M  constant.  The  value  of  this  "constant"  Is  not  known.  Try  to 
deduce  It  from  the  recorded  data  Ug(t).  You  might  have  to  make  an  Intelligent 
guess. 
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Problea  #4 


Given  arc  the  saae  facts  as  In  the  previous  problsais  Including  the 
"constant"  which  describes  the  noise  power  level. 

|fl(v)|2  «  lO**  in  -  |<  V 

Now  that  |h(v)P  is  Icnown  and  U^Ct)  Is  computable,  can  you  apply  the 
Wlcner-f liter  theory,  at  least  In  a  guessed  approximation?  Try  It  and  compute 
U4(t).  Plot  U4(t)  and  U^Ct).  Hint:  represent  |Uo(v)P  by  a  gausslan  function 
of  suitable  peak  power  and  width.  Signal  processing  specialists  always  try 
It  with  a  gausslan  function  If  they  don't  know  a  better  way. 

|0„(V)|2-P, 


(to  be  cut  by  game  monitor^ 


Problem  #5 

Try  the  same  approach  as  In  the  previous  problem,  but  with  a  guessed 
slnc^'-shaped  |Uo(v)|^. 


|Do(v)P  «  P5  slnc^(v/v5);  sine  s  - 
Plot  the  result  UsCt)  and  also  for  comparison. 


slnirt 

xs 
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Problem  #6 


Try  the  same  approach  aa  In  the  previoua  problem,  but  a  aomewhat 
different  gueaa  for  |Uo(v)P: 

-  Ps  elnc^Cv/vfi)  +  (Pq  -  P6)4oJ 

Herein  6q  meana  a  function  which  la  equal  to  1  for  v  ■  0  and  equal  to  0  for 
V  I*  0.  Plot  Ue(t)  and  Ug(t). 


(to  be  cut  by  game  monitor) 


Problem  #7 

Baa-td  on  all  of  the  accumulated  experience  try  vour  own  algnal  procaaalng 
approach  or  almply  gueaa  what  UQ(t)  might  have  been.  Call  It  U7(t).  Plot 
U7(t)  and  Uj^Ct). 
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§2.4  About  the  Signal  Synthesis 


The  noise  N(t)  covers  a  time  range  t  ■  0,  1,  2... 63  and  an  aaplltude 
range  0  <  N(t)  <  99  -  N.  The  values  for  N(t)  are  plcke'  from  a  table  of 
random  numbers.  Hence  ve  may  assume 


N(t)  -  N/2  M  50; 

N  -  99; 

N(t)  -  N  +  n(t); 

-N/2  <  N(t)  <  +N/2. 

Since  the  time  resolution  step  Is  6t  ■  1  the  frequency  range  Av  Is  Av  ■ 
l/6t  ■  1; 

The  time  duration  At  Is  At  >  64.  Hence  the  frequency  resolution  Is 
6v  ■  1/At  -  1/64.  In  other  mrds  the  Fourier  transformation  by  means  of  a 
digital  computer  provides  a  frequency  spectrum  at  v  ■  -32/64,  -31/64, .. .-1/64, 
0,  +1/64,... +31/64. 


Since  n(t)  has  a  zero-mean,  and  since  n(t)  consists  of  (almost)  unrestrained 

random  numbers,  we  may  expect  a  "white"  noise  power  spectrum. 

(|fl(v)|2)  ■  constant  •  rect(v); 

(l  ln-i<v<  +  -^ 
rect(v)  ' 


elsewhere 

The  value  of  the  constant  can  be  computed  on  the  basis  of  the  Plancherel- 
theorem: 


yifl(v)|2  dv  ■  yin(t)|2  dt. 


The  right-hand  side  will  be  computed  based  on  the  histogram-method.  It  Is 


obviously 


7’ 


|n(t) |2  dt  -  T  |n(t) |2 


The  probability  for  n  Is  constant  between  -N/2  and  +N/2,  the  constant  being 
defined  by  the  normalization  Integral 
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p(n)  ■  cooitant  *  rcct(ii/N); 


1  ■  p(n)  dn  -  N  •  conetant;  constant  ■  1/N. 

Now  wa  can  conpute  tha  aaan  aquara  of  n(t)  by  Maana  of  the  hlstograa 
nethod: 

p(n)  dn  -  (1/N)  /  n^  dn  -  N^/ll. 


|n(t)p 


Since  It  la 


■/” 


/+N/2 
•N/2 


and 


L 

/' 


|il(v)  P  dv  -  |ftP 


|n(t)  p  dt  -  T  |nP  -  T  nV12 
we  conclude  on  the  baala  of  the  Plencherel  theoraa 

|fiP  -  T  n2/12;  |flP  -  (16/3)  lO**. 

In  our  case  it  la  T  ■  64  and  N  m  lOO. 

Now  we  will  devise  the  object.  It  shall  consist  of  several  square  boxes, 
all  with  the  sasie  width  B,  hut  with  different  aaplltui'.es  An*  and  with  center 
locations  t^: 

M 

Uo(t)  -  ^  recti (t  -  t,)/BJ. 

The  number  M  Indicates  the  total  number  of  square-boxes.  He  chose  so  that 
less  than  half  of  the  time-domain  Is  covered  by  square-boxes.  The  relative 
coverage  Is  given  by: 

"relative  coverage"  ■  MB/T. 

He  will  select  the  locations  t^  such  that  there  Is  no  overlap  among  square- 
boxes.  The  object  spectrum  Is: 


Uo(v)  -  /Uo(t)  dt 


■  ^  A|,  B  slnc(v5)  e  ^**''^*. 

The  peak  power  Is: 

|Uo(0)P  - 
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At  other  frequencies  the  object  power  spectnm  Is: 

|e„(v)|2  .  b2  .Inc^CvB)  12;*, 

Thu  modulus-square  expression  can  be  re-wrltten  In  a  font  suitable  for  cosi- 
putlng  the  ex7ectatlon  value: 

-  A^j 

The  ^^(mf^J)  will  not  be  exactly  zero,  but  can  be  expected  to  be  zero  In  average 
over  the  v-domaln  (except  for  v  ■  0).  Hence  we  get  as  expected  value 
<|Uo(v)|2>  -  b2  alnc2(vB)X  V- 

He  may  also  compute  the  total  power  spectrum,  most  easily  on  the  basis  of  the 
Plancherel  theorem: 


J|Uo(v)|2  dv  -  dt  - 


RequlrCTents 

The  relative  coverage  of  the  time  domain  by  square-box  signals  shall  be 
less  than  50Z: 

MB/T  < 

The  slnc-square  function  of  the  object  power  spectrum  shall  cover  not  more 
than  half  of  the  frequency  range.  Specifically  we  request  that  the  second  zero 
of  slnc-square  Is  halfway  out  to  the  bandllmlt  1/2: 
vB  -  2;  V  -  1/4  B  -  8. 

Let  the  ratio  of  object  power  spectnim  to  noise  power  spectrum  be  about 
2  at  the  second  side-maximum  of  the  slnc-square  function.  The  second  side 
maximum  occurs  at  the  argument  5/2.  The  slnc-square  function  Is  there 
slnc2(5/2)  ■  (5ti/2)“2  «  16/1000.  Hence  we  request: 

|Uo(5/2B)|2  «  b2(16/1000)  2]  "  2  |n|2. 
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It  was  |flP  >  T  n2/2:  T  >  64;  N  >  IGO;  B  8.  Hence  we  get: 

2] A.’'  «  10  n2  -  105. 

The  relative  coverage  requirement  yields  with  B  >  8  and  T  64  for  the 
number  M  of  square-box  signals  M  <  4.  He  will  take  M  3,  which  corresponds 
to  a  relative  coverage  of  3/8. 

Now  we  select  the  object  amplitudes  Aj,  A2»  A3.  Let  us  Introduce  for 
convenience  the  term  "relative  amplitudes"  by  dividing  the  amplitudes  by  the 
mean  noise  value  N/2:  >  K^/(V/2)  >  2Ab/M.  This  reduces  the  former  soli¬ 

tude  requirement  to: 

2](N/2)2  aa^  10  N^;  40. 

He  select  now 

a^  **  6f  a^  *  2^  a^  *  1. 

It  Is  to  be  expected  that  the  square-box  #1  exceeds  the  noise  clearly »  the 
square-box  #2  Is  barely  bigger  than  the  noise »  while  the  third  square-box 
vill  be  very  difficult  to  detect*  If  at  all.  Maybe  It  Is  not  quite  as  based 
since  the  root-mean-square  noise  J(N(t)  -  N)2  •  J  |n(t)  |2  -  1^2/12  Is  less 
than  the  mean  noise  N  >  N/2. 

The  square-box  positions  t,  were  chosen  as  t^  ■■  44.5;  t2  ■■  13.5;  ts 
55.5.  The  true  original  signal  Uo(t)  Is  plotted  below. 
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§3  The  Talbot  Interferoneter 


testing  optical  components  with  the  Talbot  Interferometer  setup  resulted 
from  the  necessity  of  having  a  good  collimating  lens  for  our  scanning  Talbot 
Spectrometer.  Thus  we  have  a  sensitive  and  Inexpensive  means  for  testing  these 
components.  The  Ronchl  rulings  used  with  the  setup  give  us  the  lateral  deriva¬ 
tive  of  the  wavefront  deformed  by  the  test  object.  Thus  the  deviations  from  the 
expected  surface  shape  can  be  observed.  The  following  is  a  detailed  report  that 
explains  this  setup. 


S3.1  Introduction  to  the  Application  of  the  Talbot  Effect 

for  Measuring  Phase  Objects. 

A  phase  object  changes  the  amplitude  of  a  wave  passing  through  It  by  a 
constant  factor  but  alters  the  phase  from  point  to  point,  which  yields  a  trans¬ 
mittance  function  t(x,  y)  ■  a  e'^^^>y).  Examples  of  phase  objects  are  trans¬ 
parent  biological  objects  whose  Index  of  refraction  Is  nearly  equal  to  that  of 
Its  surrounding  medium,  and  objects  made  of  a  transparent  medium  such  as  lenses 
and  prisms.  The  techniques  now  available  to  examine  phase  objects  Include 
Zernlke's  phase  contrast  microscopy.  Interferometric  methods,  Ronchl  grating 
techniques,  Schlleren  methods  and  molrd  fringes. 

The  "Talbot  effect"^  Is  also  known  as  "Fourler-lmaglng"  or  "self-imaging". 
As  Talbot  observed  In  1836,  an  isMge  of  a  grating  Is  formed  at  a  distance  2d^/l 
(where  d  >  grating  period)  behind  the  grating,  which  has  been  Illuminated  by  a 
monochromatic  plane  wave.  The  space  between  the  grating  and  Its  Talbot  Image 
Is  empty.  In  other  words  the  Talbot  effect  Is  a  method  for  forming  Images  of 
a  periodic  object  without  any  lenses  or  mirrors.  Vfhen  properly  modified  the 
Talbot  effect  affords  a  unique  method  for  quantitatively  observing  phase  objects 
through  the  use  of  two  Ronchl  rulings,  one  of  which  Is  self-imaged  onto  the 
second.  The  complete  theory  requires  scalar  wave  diffraction  theory  and  Is 
presented  In  S3.3. 

An  explanation  using  rays  will  acquaint  the  reader  with  the  basic  opera¬ 
tion  of  the  Instrument.  However,  this  Is  only  a  crude  explanation  since  It 
does  not  take  Into  account  the  wave  nature  of  light  which  results  In  some 
possibly  useful  color  effects.  Therefore  In  the  third  section  wavelength  con¬ 
siderations  will  be  added  to  show  the  limitations  and  uniqueness  of  the  Instru¬ 
ment. 


§3.2  Qualitative  Explanation 


A  Ronchl  ruling  Is  Illuminated  with  collimated  light  as  shown  In  Figure 
2.  The  diffracted  waves  are  modulated  and  form"Fresnel  Images''^.  At  Integer 
multiples  of  the  distance  zj  2d^/X,  where  d  Is  the  grating  period,  an  Image 
of  th?  Ronchl  ruling  Is  formed.  If  a  second  grating  Is  placed  at  one  of  these 
Image  planes  In  antl-posltlon  to  the  Image,  no  light  will  be  transmitted  beyond 
the  second  grating.  Herce  a  ray  that  passes  through  a  slit  of  the  first  grating 
Is  blocked  by  the  second  grating  as  shown  In  Figure  2.  A  prism  of  wedge  angle 
a  Is  Introduced  at  a  distance  z  from  the  secb'>d  grating  and  normal  to  the 
Incident  rays.  The  rays  passing  through  the  prism  will  be  bent  through  an 
angle  c .  Thus  the  Intensity  of  the  light  as  seen  by  the  observer  varies  with 
the  angle  c  as  sketched  In  Figure  3,  with  peak  Intensities  given  by 
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e  «  tan  “  (n  +  1/2)  d/z 
and  dark  fringes  by 

e  tan  ■  md/z.  (1) 

The  angle  e  is  related  to  the  wedge  angle  a  under  small  angle  approximation  by 

e  -  (n  -  1)  a.  (2) 

A  prism's  thickness  Increrses  linearly  with  the  distance  from  its  apex,  and 
the  wedge  angle  a  giving 

a  tan  a  >  At/Ax,  where  At  *  prism  thickness 

Ax  ••  distance  from  apex. 

Now  a  phase  object  cafi  be  thought  to  be  made  up  of  prisms,  a  fact  which 
is  often  used  to  explain  the  action  of  a  lens,  so  that  for  any  point  (x,  y)  on 
the  object,  the  prism  angle  a  is  given  by 

a(x,  y)  « 

which  Ir  turn  is  related  to  e  by  substitution  into  (2) 

e(x,  y)  -  (n  -  1)  •  (3) 

The  y  derivative  3t(x,  y)/3y  is  of  no  concern  here  since  the  grating  bars  are 
parallel  to  the  y-axis. 


Thus  a  dark  fringe  structure  is  obtained  whenever  e(x,  y)  meets  the  con¬ 
ditions  of  equatl)n  (1).  Equations  (1)  and  (3)  are  plotted  in  Figure  A  as  a 
function  of  x.  Dark  fringes  occur  at  the  intersections  of  equations  (1)  and  (3). 
In  other  words  the  fringes  indicate  lines  of  equal  3t/3x,  which  is  the  x-compo- 
nent  of  the  surface  gradient  At.  Where  the  lines  are  close  together,  the  slope 
3t/3x  changes  rapidly.  Hence  the  fringe  density  or  fringe  frequence  is  propor¬ 
tional  to  the  x-component  of  the  curvature  3^t/3x^.  The  maximum  detectable 
curvature  is  determined  by  the  minimum  detectable  fringe  separation.  In  a 
sense  what  is  observed  is  the  second  derivative  of  the  phase  object. 


There  are  two  problems  of  fringe  detection.  The  first  occurs  when  the 
fringes  beoome  too  fine  for  the  eye  to  resolve  them.  This  is  only  a  practical 
limitation  and  a  magnifying  glass  would  extend  the  sensitivity.  But  ehe  secomd 
problem  is  more  fundamental.  The  field  of  view  is  finite,  say  B  ■  Nd,  where  N 
is  the  number  of  lines  in  the  grating.  For  an  object  whose  curvature  is  at  the 
minimum  detectable  sensitivity  its  curvature  must  be  large  enough  for  at  least 
one  full  (or  maybe  one-half)  fringe  across  the  field.  Hence 

[(3t/3x)„^^  -  (3t/3x)„i^J  (n  -  1)  >  d/2z. 

If  (3t/3x)^^y  is  at  X  “  +  Nd/2  and  (9t/3x)jjjj^jj  is  at  x  -  -Nd/2,  then 

3t(x  +  Nd/2,  y)/3x  -  3t(x  -  Nd/2,  y)/3x  «  32t/3x2  Nd. 


Hence, 


|d2t/3x2|  > 


d/2 

z(n  -  1)  Nd 


1 

2Nz(n  -  1)  • 


(4) 


The  maximum  Talbot  distance  is  given  by 
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-  Hd2/2X. 

Using  this  as  z  In  equation  (4)  yields 

|32t/3x2|  >  -  1)  '  B^^n  -  1) 

For  exaaple  If  X  ■  5  x  10~'*  an,  n  -  1  ■  1/2,  B  ■  10^  an  we  obtain  l32t/3x2|  > 

10-7  .  1/10  ksi. 

As  an  exaaple  of  the  system  a  perfect  thin  lens  Is  considered.  The  phase 
of  this  lens  Is  given  by^ 

<P(x.  y)  -  (-ir/XfXx2  +  y2) 

where  f  Is  the  focal  length  of  the  lens,  and  (x,  y)  Its  coordlnatas.  The 
relative  phaae  of  any  object  la 

9(*.  y)  "Qit/XXn  -  1)  t(x,  y) 

.  .  t(x,  y)  -  -(x2  +  y2y(2f(n  -  1)1, 
and 

3t(x,  y)/3x  -  -  x/Kn  -  l)f) 
which  when  substituted  Into  equation  (3)  gives 
e(x,  y)  -  -x/f. 

This  Is  a  straight  line  which  whan  plotted  onto  Figure  4  gives  equally-spaced 
fringes  of  fd/z.  An  Imperfect  lens  will  not  have  straight  equally  spaced  fringes* 
whereby  Its  quality  can  be  dAtarmlned. 

It  Is  worth  noting  that  this  system  likewise  tests  ths  collimator  objective 
and  can  serve  as  a  good  method  for  testing  lenses.  In  addition  the  ability  to 
completely  darken  the  field  behind  the  second  grating  Is  a  test  of  the  degree 
of  colllMtlon  of  the  objective  since  converging  or  diverging  waves  produce 
aralri  fringes  behind  the  second  grating  due  to  unequal  grating  periods. 

The  work  of  Nlshljima  and  Oster  can  be  explained  by  rays?  However,  on 
the  other  hand  a  complete  explanation  of  our  Instnaent  requires  the  use  of 
diffraction  theory  which  accounts  for  the  color  effect  and  which  explains  the 
ultimate  limitation  of  this  Instrument.  The  location  of  the  grating  Image  Is 
wavelength  dependent,  a  fact  that  has  been  exploited  In  designing  filters  and 
a  spectrometer^.  Furthermore,  the  complete  setup  (vide  Figure  5)  allows  the 
selection  of  the  first  or  second  derivative  of  the  phase  object. 
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Figure  4.  Plot  of  equaclone  (1)  and  (3).  Equation  (1)  la  the  aet  of  hori¬ 
zontal  !<tralght  lines.  X  aarka  the  locations  of  the  fringes  for 
this  test  object. 
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Figure  5.  Test  setup  Incorporating  a  spatial  filter  plane  for  selecting 
first  or  second  derivative. 
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§3.3  Theoretical  Explanation 


The  ray  explanation,  which  la  eaaentially  due  to  Nishijlaa  and  Oster  and 
is  presented  in  §3.2,  helped  one  to  understand  the  basic  operation  of  the  instru- 
nent.  In  this  section  a  rigorous  derivation  is  presented  based  on  scalar 
diffraction  theory.  The  result  will  be  an  interpretation  of  the  Inage  as  two, 
three  or  many  shifted  object  wave  fronts,  soaetiaes  tilted.  Accordingly  we 
will  use  the  teras  "shearing  interferoaetry"  as  is  coaaon  when  two  shifted 
object  wavefronts  interact.  In  generalization  thereof  we  will  introduce  the 
teras  "triple  shearing  Interferoaetry"  and  "aultiple  shearing  interferoaetry". 
In  soac  instances  these  shearing  interferences  will  represent  the  first  or  the 
second  derivative  of  the  object.  The  instruaent  is  diagraaaed  in  Figure  6, 
which  illustrates  the  pertinent  distances. 


Consider  a  aonochroaatic  plane  wave  incident  on  grating  G1  in  the  plane 
at  -z^.  The  field  behind  this  grating  is 


u(x,  y,  -Zj) 


-Ikzi  r  -  2ninx/d 

e  Z  Cn  « 

n"-“ 


where  the  periodic  grating  is  considered  infinite  in  extent  and  of  period  d, 
and  is  expressed  as  a  Fourier  exponential  series;  k  ■  2ir/X  where  X  ■  wavelength 
of  light.  Following  a  line  of  reasoning  like  Edgar's^,  the  incident  plane 
wave  is  by  action  of  the  grating  replaced  by  a  set  of  plane  wavefronts  whose 
x-direction  cosines  are  nl/d  and  whose  field  strength  is  proportional  to  C^. 
This  gives  a  propagating  field  behind  G1  as 


,  .  'Ikzi 

u(x,  y,  z)  «  ft  ^ 

which  at  the  object  plane  z  ■  0 

u(x,  y,  0-)  - 


r  ^  ^lk[nXx/d+(zi+z)/l-(nX/d)^) 

^  ''n 

n 

reduces  to 

^  ^  ^ikCnAx/d+zi/l-inX/d)^] 


n 


Applying  the  Kirchoff  boundary  conditions  to  the  field  as  it  passes 
through  the  object  as  we  did  through  grating  Gl,  we  have 

/  -ikz;  V  n  t  \  ik  [  nXx/d+z  j  /I  -  (nA  /d )  ^  ]  ... 

u(x,  y,  0+)  -  e  ^  i  Cp  uo(x,  y)  e  ‘  *  (5) 

n 

where  uqCx,  y)  is  the  two  dlaensional  object  transaittance  function. 

Each  of  the  plane  waves  is  diffracted  as  it  passes  through  the  object. 

At  this  point,  the  concept  of  the  angular  spectrum  is  very  useful.^  In  this 
formulation  the  diffraction  phenomenon  is  a  multiplicative  quadratic  phase 
factor  exp(ikz/l  -  12 (v^  +  pi)  )  in  the  Fourier  domain  which  increases  with  the 
propagation  distance  z.  Let  0(v,  p;  (H-)  be  the  angular  spectrum  of  equation 
(5)  defined  by 

a(v,  u;  0+)  -  |j  u(x,  y,  0+)  e~2iri(vx+uy) 

where  the  integration  in  this  section  will  be  over  the  Interval  (-»,  <•)  when 
unspecified,  and  (v,  u)  are  the  Fourier  spatial  frequency  components.  The 
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angular  spectrum  for  the  plane  Just  preceding  G2  (z  -  Z2-) 


0(v,  p;  zj-) 


uCv,  w;  CH-) 


Ikz  2 /l-X^ (v^+w^) 


Is 


-  I  C„  <I„(V  .  n/d,  d)  .lMr.A-(„i/d)2«2/l-X^(v24v2)) 

n 

where  OqCv,  p)  Is  the  angular  spectrum  of  the  object.  The  Klrchoff  boundary 
conditions  are  applied  to  the  second  grating  at  z  ■  Z2. 


u(x,  y,  Z2+)  -  u(x,  y,  Z2“)  I  e^™  g2iilmx/d 
Itto  ^ 

where  the  e  factor  accounts  for  the  fact  that  the  second  grating  la  shifted 
by  half  a  period  with  respect  to  Gl.  In  the  Fourier  domain  this  becomes  a  con¬ 
volution  which  results  In 


a(v,  p;  Z2+)  -  e  I  I  “o(''  ~  (n^)/d.  y) 

”  *  ^lk{zi A-(nX/d)^  +  Z2 A-X^[ (v-m/d)^+p^]} 

(6) 

The  Image  plane  Is  conjugate  to  the  object  plane  z  >  0.  Therefore  we 
compute  the  field  virtually  back  to  the  z  ■  0  plane  In  order  to  find  the 
resultant  field  In  the  plane  of  observation,  and  give  the  field  a  new  symbol 
v(X|  y)  with  A^(v,  p)  as  Its  Fourier  spectrum 

_/  s  .\  “lkz2/l”X2  (v2+p2) 

^(v,  p)  ■  Q(v,  p;  Z2+)  e  ^ 


Substituting  equation  (6)  Into  this  and  simplifying  by  using  the  first  two 
terms  In  the  Taylor  serlas  expansion  of  the  square  roots  In  the  exponential 
terms  yields 

«(v.  li)  -  n  Cn  c.  ao(''  -  (■rtO/'l.  a)  -Xl<n»/a)2-l2<11ix/d)2l_ 

n  m 

But  as  explained  In  the  first  section  the  two  gratings  are  separated  by  some 
multiple  M  of  the  Talbot  distance  zj.  Thus  putting  Z2  +  zi  ■  M  Zj  m  2Md^/X 
Into  the  above  equation  reduces  It  to 

^(v,  p)  -  I  I  Cn  Cjj  e^™  Qo(v  -  (n+m)/d,  p) 

”  ®  ^2xl[Xz2vm/d-Z2X(m^-n2)/2d^] 

Taking  the  Fourier  transform  of  this  gives  the  resultant  field  In  the  Image 
plane. 

V<x,  ,)  ■  n  C„  C.  u„(x  +  .xjX/d,  y)  e2>ll»/2+(.«)2x2/xT*(M.)x/dJ 
n  m  '  (8) 


Thus  multiple-shearing  Interferences  are  observed  In  the  linage  plane.  In 
ordinary  shearing  Interference  only  two  terms  are  present,  u(x  -f  Ax,  y)  -  u(x,  y) . 
In  order  to  reduce  the  large  number  of  terms  of  the  double  series  we  will  perform 
certain  spatial  filtering  operations.  The  field  expressed  by  equation  (7)  Is 
found  In  the  spatial  filter  plane  located  one  focal  length  behind  the  first 
lens  of  the  telecentrlc  system.  In  order  to  filter  Just  the  n  -f  m  ■  0  or 
n  -f  m  -  1  term,  which  we  will  show  leads  to  special  Interferences,  we  must 
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Figure  6.  Talbot  Interf eroneter  setup.  Telecentrlc  systen  is  focused  onto 
the  object.  The  spatial  filter  plane  It  used  to  slapllfy  the 
Bultlple-shearing  Interferences. 


Figure  7.  The  spatial  filter  plane  with  an  arbitrary  object  spectrum.  The 
spectrum  Is  shifted  along  the  v-axls  by  1/d. 
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assume  that  the  shifted  spectrums  ~  (n^)/d,  p)  do  not  overlap.  This  is 

shown  In  Figure  7  along  with  the  bandwidth  restriction  of  the  object,  |Av|  5 
l/2d,  l.e.,  i.q(v,  p)  «  0  whenever  |v|  ^  l/2d. 


§3.4  Zero  Order  Filtering 


The  multiple  shearing  Interferences  can  be  simplified  by  spatial  filtering. 
Our  first  consideration  Is  the  zero  order  filtering  case,  n  -f  m  ■  0.  Thus 
equation (8)  reduces  to 


v(x,  y)  -  I  C  C  e  ^  uo(x  +  mzX/d,  y)  (9) 

m 

where  Z2  Is  written  as  z  without  the  subscript  2  and  Is  the  distance  from  the 
object  to  the  second  grating.  If  we  restrict  ourselves  to  only  considering 
either  amplitude  or  phase  modulated  gratings  that  are  symmetrical  about  the 
origin  then  Cg,'s  are  real  and  *  C-q.  Equation  (9)  becomes 

v(x,  y)  -  Co^{uo(x,  y)  +  I  (Cbj/Cq)^  uqCx  +  mzA/d,  y)l. 

^  mftO  } 

These  multiple-shearing  Interferences  can  be  simplified  to  triple-shearing 
Interferences  If  (Cq/Cq)^  terms  for  all  but  the  first  pair  are  negligibly 
small.  By  using  a  special  grating  whose  transmittance  Is  given  by  g(x)  «  1  -f 
cos(2iTx/d),  triple-shearing  Interferences  would  result,  since  Cq  "  0,  |m|  >  2. 
Alternatively,  a  grating  can  be  designed  which  is  still  binary  and  yet  meet 
this  condition.^  We  now  consider  the  use  of  an  inexpensive  Ronchi  ruling  as 
the  grating.  The  Fourier  coefficients  are  Cg,  ■  1/2  slnc(m/2);  all  even  Cq  ■  0, 
where  slnc(x)  ■  sln(iTx)/itx.  Thus 


v(x,  y)  -  Co^|uo(x,  y)  I  slnc2(m/2)  uqCx  +  mzA/d,  y)|.  (10) 

(odd) 

These  multiple-shearing  Interferences  can  be  approximated  by  triple-shearing 
Interferences  In  this  case  since  the  |m|  ■  3,  5,  7,...  terms  can  be  neglected. 
The |m|  -  3  term  Is  only  one-ninth  that  of  the  |m| ■  1  term.  Higher  order  terms 
are  much  smaller  than  this.  Therefore  equation  (10)  reduces  approxlnately  to 

v(x,  y)  l/4|uo(x,  y)  -  (2/v)^  uo(x  +  zA/d,  y) 

-  (2/ii)2  uo(x  -  zA/d,  y)|.  (11) 


These  triple-shearing  Interferences  will  sometimes  be  roughly  the  second 
derivative,  3^u/dx^  of  the  object  If  the  shift  zA/d  Is  smaller  than  the  resolu¬ 
tion  length  6x,  which  Is  the  Inverse  of  the  object  bandwidth  Av  ■  l/6x. 

Actually  zA/d  ■  (l/3)6x  Is  still  good,  as  we  will  show  by  making  use  of  the 
Bernstein  theorem.  First  we  expand  the  last  two  terms  In  equation  (11)  In  a 
Taylor  series  around  zA/d. 

00 

u(x  +  zA/d,  y)  •  I  Uq*’^(x,  y)  1/p!  (zA/d)*^ 
p-O 


Uq(x  -  zA/d,  y)  ~  I  Uo*’^(x,  y)  (-l)*^/pl  (zA/d)*^ 
p-O 
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i 


where  uq^^x,  y)  ■  d^uo(x,  y)/3x^.  Therefore  equation  (11)  becoaes 

f  * 

v(x,  y)  -  l/4|uo(x,  y)  -  (2/it)2  I2uo(x,  y)  -  2  I  uo**\x,  y) 

(  p-2 

1/pl  (zX/d)Pl|. 

The  suomatlon  term  can  be  approximated  by  the  second  derivative,  l.e., 

00 

[  uJP^x.  y)  1/pl  (zA/d)P«  up^(x,  y)  l/2(zA/d)2  (12) 

p-2 

(even) 


if  the  terms  p  -  4,  6,..  are  small  compared  to  p  -  2.  This  can  be  shotm  as 
follows.  Let 

fp  -  1/pl  (zA/d)**  uo**\x,  y). 

Then  equation  (12)  becomes 


If 


I 

p-2 

(even) 


f 

P 


M  fj 


ifpi  ;si/io  fr*. 

The  factor  of  1/10  Is  picked  as  a  figure  of  merit,  but  other  values  may 
be  used.  As  we  shall  see,  this  factor  determines  the  maximum  distance  that  the 
object  can  be  placed  from  the  second  grating  G2.  Now, 

|f2l  2  f?**  -  l/2(zA/d)2  (13) 

By  the  theorem  of  Bernsteln^^,  the  derivatives  of  a  bandllmlted  function  are 
bounded  by  the  maximum  value  that  the  function  attains.  That  Is, 

|g^P\x)|  <  (2irL)P  Bo 

where  Bq  -  !(8(x))aiax  ^  maximum  frequency  of  g(x)  In  the  Fourier 

domain.  Call  2L  -  Av;  we  have 

|f  I  <  1/pl  (zA/d)  (irAv)***^  ■  2/p  I  (zAirAv/d)**”^  f“* 

p  bbx 

which  when  compared  with  equation  (13)  Implies  that 
2/pl  (zAirAv/d)**"^  <  1/10, 
zA/d  <  (pI/20)^^P'^/(irAv). 

The  factorial  dominates  the  root  process  as  can  be  shown  using  the  Stirling's 
approximation  for  factorial, 

(pI/20)^^**'”^  M  p/®  large  p. 
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Therefore  Che  smallest  p,  (p  ■  4)  is  Che  cightesc  bound  for  z,  and  we  have 

zX/d  ;s  /2472O  /  irAv  «  6x/3  .  (14) 

Recall  chat  for  spatial  filtering  we  require  that  4v  <  1/d.  Therefore  the  most 
stringent  condition  is  for  4v  -  1/d.  This  is  required  in  order  that  the  shifted 
object  spectrum  will  not  overlap, as  shown  in  Fi' are  6.  Thus  the  object  to  G2 
distance  z  in  this  instance  must  be  less  than  about  1/6  of  a  Talbot  distance, 
or  more  exactly, 

z  <  1.1  z^/2ii.  (15) 

Under  the  above  restrictions  as  given  by  either  equation  (14)  of  (IS)  for  the 
maximum  distance  z,  the  field  in  the  image  plane  is  given  by 

v(x,  y)  «  1/4(1  -  S/ir^)  {uo(x,  y)  -  _  2  u^^x,  y)} 

«  1/20  {uo(x,  y)  +  5/7  (zX/d)^  up^(x,  y)}.  (16) 


§3.5  First  Order  Filtering 


The  mathematics  for  this  case  follows  simply  from  the  previous  section. 
Here  we  filter  either  the  n+B"lorn'fm*-l  term.  The  only  difference  is 
in  the  sign  of  a  phase  factor  which  is  unimportant  when  detecting  the  signal. 

The  field  in  the  image  plane  becomes  for  the  n  -f  m  ■  +1  case, 

v(x,  y)  -  Co  Cj  {uo(x,  y)  -  uo(x  +  zX/d,  y)}  (17) 

This  surprising  result  is  obtained  since  the  even  coefficients  are  zero.  Those 
terms  which  contribute  under  the  n  -t-  m  -  1  condition  are  Cg,  Ci,  ^  Cg  Cg|.i.  This 
combination  implies  that  for  any  m,one  coefficient  will  always  be  even  and  henee 
equal  to  zero.  This  is  very  much  like  ordinary  shearing  interferometry  and  we 
will  again  show  that  under  a  certain  condition  this  is  approximately  the  deriva¬ 
tive  of  the  object.  We  ignore  the  phase  since  it  is  the  intensity  which  is 
observed  in  the  image  plane,  and  deal  with  the  bracketed  terms.  The  last  term 
is  expanded  in  a  Taylor  series  around  zX/d  and  when  combined  with  uo(x,  y)  gives 

I  u^^’^x,  y)  (zX/d)**  1/pl  wuo^^(x,  y)  (zX/d) . 

Pil 

Using  the  same  kind  of  arguments  as  in  the  previous  section  this  approximation 
is  valid  if  the  shift  is  less  than  1/15  5x  (equation  18). 

zX/d  1/5ttAv  «(1/15)  6x.  (18) 

Under  the  strictest  condition  |4v|  5  1/d  the  condition  of  equation  (18)  means 
that  the  object  to  G2  distance  is 

z  S  zf/lOn.  (19) 

Therefore  the  field  in  the  image  plane  for  first  order  filtering  under  the 
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conditions  of  equation  (18)  or  (19),  whichever  may  apply,  Is  given  by 

v(x,y)  l/ir  (^^/d)  y) .  (20) 

One  other  case  has  been  under  Investigation.  Here  we  record  the  Intensity 
pattern  on  film  or  on  a  diffuser  (like  a  rear-project Ion  screen),  then  view  It 
with  a  low  pass  detector  such  as  the  eye.  The  experimental  results  were  shown 
In  the  previous  case  where  no  filtering  was  done.  The  reader's  eye  performed 
the  filtering  as  he  observed  the  results.  Unfortunately  the  mathematics  has 
proven  cumbersome.  But  we  hope  to  find  the  necessary  approxlaiatlon  to  yield 
a  simple  Interpretablc  solution. 


§3.6  Experiments 


The  Talbot  Interferometer  for  testing  phase  objects  was  placed  on  a 
laboratory  optical  bench.  It  used  two  Ranchi  rulings  of  10  lines  per  mm. 

Although  they  were  quite  scratched,  they  performed  satisfactorily  because  the 
resultant  Intensity  In  the  plane  of  observation  Is  summed  over  a  large  number 
of  lines.  One  grating  was  mounted  for  rotational  aiotlon  while  the  other  had 
translational  movement.  These  two  motions  allowed  the  gratings  to  be  aligned 
as  desired.  The  gratings  were  held  in  bench  carriers  which  accomplished  longi¬ 
tudinal  motion. 

Both  He-neon  laser  and  white  light  sources  were  used.  The  laser  beam  was 
focused  onto  a  pinhole  and  thereafter  collimated  (diameter  SOm) .  For  the 
white  light  source  exper;lments,  a  500  wa'.t  slide  projector  was  used.  To  Increase 
the  throughput  an  adjustable  slit  was  ur.ed  Instead  of  a  pinhole;  It  was  aligned 
parallel  to  the  bars  of  Gl.  The  width  of  the  silt  was  then  made  sufficiently 
narrow  so  that  the  self-images  trould  appear  with  good  color  contrast.  It  was 
discovered  that  the  molrd  fringes  obtained  behind  the  second  grating  when  both 
gratings  were  aligned  was  a  measure  of  the  degree  of  colliiutlon.  This  work 
has  been  accepted  for  publication  and  the  manuscript  Is  provided  In  Appendix  B. 

The  method  of  observing  the  Interference  fringes  was  to  project  the  Image 
of  the  object  with  Its  fringe  structure  onto  a  screen.  The  phase  objects  tested 
were:  various  lenses,  prisms,  plastic  bags  and  boxes,  the  change  of  index  of 
refraction  caused  by  heat  of  a  candle  flame,  and  drops  of  plastic  resin  or  var¬ 
nish  arranged  Into  different  configurations  on  glass  slides. 

Photographs  of  these  experiments  were  made.  Two  of  these  appear  In  our 
paper  that  appeared  In  Optics  Communications,  February  1971.  A  copy  of  this 
paper  is  provided  In  Appendix  A.  The  photographs  are  of  a  candle  flame  that 
was  placed  between  gratings  separated  by  about  12  cm.  To  fully  appreciate  the 
filtering  aspect  of  this  Instrument  we  have  separated  the  gratings  by  about 
30  cm.  with  the  flame  placed  close  to  Gl.  In  this  case  the  multiple-shearing 
Interferences  are  readily  observed  by  the  great  number  of  candle  wicks  In  the 
field  of  view.  By  allowing  only  the  zeroth  diffracted  order  to  propagate  beyond 
the  spatial  filter  plane,  triple-shearing  interferences  are  obtained,  as  veri¬ 
fied  by  th'3  presence  of  the  wick  which  appears  at  three  separate  locations 
corresponding  to  the  amount  of  shear  Introduced.  Lateral  shearing  Interferences 
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are  obtained  when  first  order  filtering  Is  used,  as  Indicated  by  observing 
only  two  wick  Images.  Photographs  of  this  along  with  some  beautiful  color 
pictures  for  a  white  light  source  were  shown  at  the  1971  Annual  Spring  Heetlng 
of  the  Optical  Society  of  America  th^t  was  held  In  Tucson,  Arizona.. 


§4  The  Talbot  Autocolllsiator 


We  found  that  colllmatlon  of  the  source  was  critical  when  using  the 
Talbot  Interferometer  setup  (93).  Otherwise  It  was  not  possible  to  completely 
darken  the  field.  By  Increasing  the  separation  of  the  gratings  to  Its  SMxlmum 
permissible  extent,  a  very  sensitive  means  for  accurately  locating  a  point 
source  In  the  focal  plane  of  the  lens  was  obtained.  Thus  the  problem  of 
accurately  and  Inexpensively  collimating  the  light  sources  for  use  with  the 
Fourier  Spectrometer  was  solved.  A  manuscript  on  "A  Simple  Interferometric 
Method  of  Beam  Colllmatlon"  has  been  accepted  for  publication  In  Applied  Optics 
this  August.  The  manuscript  Is  reproduced  In  Appendix  B. 
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§5  The  Talbot  Interferometer  with  Circular  Gratings 

Earlier  (§3),  we  discussed  the  use  of  the  Talbot  Interferometer  for 
testing  optical  costponenta.  This  setup  ^<^ed  Ronchl  rulings  and  displayed  the 
derivative  of  the  phase  objci'.t  along  a  direction  transverse  to  the  ruling  axis. 
Thus  the  testing  of  an  axially  symmetric  object  such  as  a  lens  requires  that 
It  be  rotated  for  complete  examination.  Ideally,  such  objects  are  well  suited 
for  testing  with  Interferometers  that  give  a  constant  shear  to  the  wavefront 
In  the  radial  direction.  Therefore  we  intend  to  modify  the  Talbot  Interfero¬ 
meter  by  replacing  the  straight  line  gratings  with  circular  gratings.  In  this 
way  the  radial  gradient  of  the  deformed  wavefront  can  be  observed.  A  circular 
grating  (see  Fig.  9)  Is  defined  In  this  report  as  consisting  of  a  number  of 
concentric  circular  rings  of  constant  radial  spacing  on  a  transparent  surface. 
Our  description  will  be  brief  since  the  experiments  and  theory  of  this  new 
modification  are  not  quite  complete.  We  present  here  our  analysis  on  whether 
the  circular  grating  i^ll  self-image  i^en  Illuminated  by  a  plane  wave. 


§5.1  Theory  of  the  Talbot  Effect  In  Cylindrical  Coordinates 


The  theory  for  the  Talbot  set  of  self-imaging  objects  In  cartesian  coor- 
dlnites  has  been  presentrd  In  a  paper  by  U.  Montgomery. What  Montgomery 
essentially  did  was  to  solve  the  wave  equation  while  asking  the  follotrlng 
question:  "What  are  the  necessary  and  sufficient  conditions  that  the  object 
must  satisfy  In  order  that  a  faithful  Image  of  It  be  found  In  a  parallel  plane 
z  ■  d  >  0?"^2  Our  solution  uses  this  approach. 


We  want  to  solve  the  wave  equation  In  cylindrical  coordinates  for  the  set 
of  axially  symmetric  objects  of  finite  aperture  which  are  periodic  In  the 
direction  of  wave  propagation,  the  z-axis  of  Figure  8.  Therefore  we  have  the 
following  boundary  conditions: 

u(r,  <p,  z)  >  u(r,  z)  axial  symmetry  (1) 

u(r,  Az)  >  u(r,  0)  periodicity  along  z-axls  (2) 

u(r,  0)  >  0  for  r  2  1  finite  aperture  (3) 

where  u(r,  0)  Is  the  Talbot  set  of  objects  we  desire,  Az  Is  the  period  of  the 
repeating  Image  and  the  cylindrical  coordinates  (r,  9,  z)  are  related  to  the 
cartesian  coordinates  (x,  y,  z)  by 

x2  +  y2  ■  r^;  X  ■  r  cos  9;  y  ■  r  sin  9* 

The  wave  equation  In  cylindrical  coordinates  Is 


0 


where  u  -  u(r,  9  z) ,  k  ■  2ir/X,  and  X  Is  the  wavelength  of  the  light  source. 
Applying  the  axial  symmetry  condition  equ.  (1)  to  the  above  equation  reduces 


It  to 


3r2 


r 


3u  .  3^u  .  ,2 

37  3?r  ^  “ 


-  0 


(A) 
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(5) 


Next  the  separation  of  variables  technique  Is  employed,  letting 
u(r,  z)  ■  R(r)  Z(s)  ■  RZ. 

The  boundary  conditions  of  equ.  (2)  anJ  (3)  become 

R(l)  -  0  (6) 

Z(Az)  -  Z(0).  (7) 

Inserting  equ.  (5)  Into  (4)  we  obtain 

R"Z  +  -  R'Z  +  RZ"  +  k^RZ  -  0. 
r 

Rearranging  this  and  dividing  by  RZ  we  have 

^  +  -  —  +  k2  -  constant  "  -  J-  (8) 

The  left  and  right  hand  sides  of  this  equation  are  equal  to  a  constant 
since  they  are  respectively  functions  of  r  and  z  only.  If  we  let  the  constant 
be  given  by  (ky)^  and  solve  the  right  hand  side  of  equ.  (8),  we  have 

Z"  .  -  k^Y^Z 

which  yields  the  solution 
Z(z)  -  Z(0) 


The  constant  Is  chosen  positive  and  real  (k^y^  >  8)  to  avoid  evanescent 
waves  that  are  attenuated  within  a  few  wavelengths  along  the  z^axls.  These 
waves  would  not  be  present  over  distances  that  our  Instrument  will  operate. 
Thus  our  solution  Is  a  wave  propagating  In  a  direction  given  by  the  direction 
cosine,  y  (y  <  1).  Furthermore  we  choose  the  positive  sign  in  the  exponent 
since  forward  propagation  is  of  Interest  here. 


Now  the  boundary  condition  of  equ.  (7)  Is  applied  to  this  solution. 


giving  us 


lky4z 


mX  ^  , 

or  Ym  -  ^  ^  1. 


This  last  Inequality  comes  from  our  earlier  statement, 
maximum  value  that  m  can  obtain, 

1 


m 


max 


Yl  ■ 


Hence  there  Is  a 


Thus  our 


where  0  <  e  <  1  Is  a  positive  number  that  makes  nig,^  an  Integer, 
eigenvalues,  y^,  form  a  discrete  set  due  to  the  self-imaging  requirement. 
Consider  now  the  left  hand  side  of  equ.  (8)  which  when  multiplied  by  R  yields 


I"  +  7  R'  +  k^d  -  y2)r  . 


0. 


Let  C  ■  r  k  /l  -  y^  and  substitute  It  Into  the  above  equation.  We  have 

dlMil  +  1  Mil  +  R(r)  -  0 

This  Is  the  Bessel's  differential  equation  of  order  zero  that  has  the  solution 
R(C)  -  Jo(kr  /l“  -  72) 

Applying  the  boundary  conditions  of  equ.  (6)  to  this  solution  gives 
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R(l)  -  0  -  JqCv,) 


where  Vg|  ■  k  A  -  y2  are  the  zeros  of  Jg.  Thus  our  eleaentary  solution  Is, 
to  within  a  constanf,  given  by 

u,(r,  z)  -  Jo(v*r)* 


Our  total  solution  Is  a  superposition  of  all  our  separate  eigenvalue 
solutions.  ■max 

u(r,  *)  -  I  C,  jQ(kr/l  -  (nYi)^  )  (9) 

Bt*0 

In  order  to  evaluate  the  Cb'*  use 


(10) 


where  Vq,  Vq  are  zeros  of  Jq,  and  6^  Is  the  Kronecker  delta, 

.  .  /  0,  n  ^  a 

na  \  1 ,  n  ■  a 

Our  Initial  condition  Is  u(r,  0),  that  Is  u(r,  0)  Is  rhe  Talbot  sat  of  objects 
with  axial  syaaetry.  Proa  equ.  (9) 

u(r,  0)  -  J  C,  Joe'll).  (11) 

a-0 

By  Integrating  this  as  follows 


i: 


u(r,  0)  Jo(vnr)  rd* 


the  Cgi'a  can  be  deteralned.  Substituting  equ.  (11)  Into  the  above,  and 
exchanging  susautlon  and  Integration  yields, 

[  u(r,  0)  Jo(Vnr)  dr  -  (^  [  Jo(%r)  Jo(Vnr)  rdr 

)o  a-0  Jn 


■0 

-  I  C,i6^Jl2(%) 

a-0 


0 

2  ^na 
T  ^n  • 


Therefore, 


2 

“  Ji'2(V  J  rdr. 


Thus  our  solution  for  the  Talbot  set  of  axial  syaaetrlc  objects  Is: 

*aax  rl 

u(r,  *)  -  I  C„  Jq(>^),  Cb  -  -j-  TT  0)  Jo('^*^)  'dr 

B-0  ^  ^0 

where  a,..,  -  ^  -  e.  Vb  "  k/1  -  (byi)2. 
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SS.2  Self-Imaging  of  a  Circular  Grating 

The  circular  grating  (see  Fig.  9)  can  be  expressed  in  a  Fourier  Series  as 
follows: 

0(r)  .  I  C. 
m"0 

where  ■  a  sine  no,  a  is  the  period  and  a  is  that  fraction  of  the  period  for 
which  the  grating  Is  transparent.  Our  exact  solution  for  the  Talbot  set  of 
axially  syonetrlc  gratings  Is  In  Bessel  functions  ifhlch  for  large  arguments 
behave  approximately  like  a  cosinusoidal  function. 

JoCvn*")  “  ^  cos(Vnr  -  >  25.  (12) 

Therefore  our  elementary  solution  In  r  will  build  up  to  a  periodic  structure 
like  a  circular  grating  if  -  a/2irR,  l.e.  k:4  -  (mYi)2  -  a/2irR. 

However,  the  l//t  factor  In  equ.  (12)  means  that  the  contrast  will  dimin¬ 
ish  with  Increasing  r.  That  Is,  an  object  which  belongs  to  a  Talbot  set  of 
axially  syametrlc  objects  must  have  decreasing  values  of  transmittance  as  r 
Increases  In  order  to  be  faithfully  self-imaged.  Thus  the  circular  grating  only 
approximately  belongs  to  this  Talbot  set  of  objects. 

Circular  gratings  were  suide  on  AGFA  SCIEMTIA  10E7S  glass  photi graphic 
plates.  These  were  used  to  replace  the ‘Ronchl  rulings  of  our  Talbot  interfero¬ 
meter.  The  Talbot  Images  observed  showed  a  decreased  contrast  toward  the 
grating  edges  which  can  be  explained  by  the  1/^  factor  mentioned  previously. 
This  may  limit  the  maximum  size  of  the  object  that  can  be  tested.  Ue  have 
looked  at  SO  mm  diameter  objects  without  difficulty.  A  full  report  can  not  be 
made  at  this  time  as  the  wor).  Is  not  quite  complete. 
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y 


Figure  8. 


The  geometry  and  coordinates  used  in  §5. 
and  has  values  u(r,  0). 


Aperture  is  circular 


Figure  9. 


The  circular 
shown  on  the 


grating.  The  normalized  transmittance  function  is 
right  side. 
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The  molrd  efiact  ie  very  eeoeiUve  in  detecting  email  dlfferencee  between  two  eimilar  gratinge.  Tbeee 
dlfferenoee  might  be  caueed  by  an  object  with  pbaee  gradient  placed  between  the  two  gratinge.  The  perform¬ 
ance  of  euch  an  inetniment  can  bo  underetood  in  terme  of  the  ■Talbot  effect*  (aleo  called  ■Fourier 
imaging”  or  ■eelf-imaging').  Slight  modiflcatiooe  provide  ehearlng  Interferenoee  end  the  eecond  deriva¬ 
tive  of  the  object. 


The  moird  effect  reveals  very  small  imperfec¬ 
tions  of  two  gratings  placed  on  top  of  each  other. 
Lord  Rayleigh  used  t^  effect  for  testing  dif¬ 
fraction  gratings.  Hw  same  basic  idea  has  been 
used  for  many  other  purposes.  For  example,  in 
electron  microscopy  the  imperfections  in  two 
pieces  of  crystal  lattice  can  be  made  visible  in 
this  way  [1].  Essentially  the  same  moird  effect 
is  utilized  for  studying  the  shape  of  a  diffusely 
reflecting  surface,  when  the  s^ow  of  a  grating 
falls  onto  that  surface.  The  grating  shadow  is 
observed  through  the  same  grating  [2].  Also  the 
shape  of  a  refracting  object  can  be  investigated 
by  means  of  the  molrd  effect  [S].  The  refracting 
object  is  placed  before  or  behind  a  first  'grating 
(fig.  1).  The  shadow  of  the  first  grating  will  be 
deformed  due  to  the  refractive  gradient.  The 
moird  fringes  observed  behind  the  seconc'  grating 
placed  at  a  distance  s  from  the  object  are  lines 
of  eqiul  deviation  [eq.  (1)]. 


Fig.  1.  The  Talbot  Interferometer.  Gl,  G2  gratinge 
with  period  d\  ray  deflection  angle  C  In  object  O. 


y)  ■  wrf/*  :  m  ■  0,  ±  1 ,  ±  2. . .  .(1) 

We  intend  to  extend  this  method  of  observing 
refractive  gradients  by  means  of  moird.  As  de¬ 
scribed  so  far  this  method  is  based  entirely  on 
rays.  This  point  of  view  is  not  satisfactory  since 
it  is  known  that  wave  optical  color  effects  occur 
when  white  light  passes  through  two  gratings  at 
a  finite  distance.  This  happens  with  uncollimated 
[4]  and  with  collimated  It^  fS].  Such  color  ef¬ 
fects  may  well  prove  to  be  useful.  A  more  im¬ 
portant  (Ejection  against  the  ray-optical  point  of 
view  is  the  inability  to  explain  the  fundamental 
limitation  of  this  method.  Based  on  eq.  (1),  one 
might  suspect  that  an  arbitrarily  small  deflec¬ 
tion  angle  can  be  detected  if  only  the  distance  i 
of  the  two  gratings  is  sufficiently  large.  Since 
such  unlimited  detectibility  is  never  attainable 


Fig.  2.  Overlapping  diffraction  orders  behind  a  grating 
of  finite  width  B  ^  Nd. 
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(a) 


(b) 

Fig.  3.  Talbot  Interferencee  wltb  a  candle  flame  aa  object  between  the  two  gratings,  a.  Gratings  parallel,  b.  One 

grating  slightly  rotated. 
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one  might  try  to  patch  on  the  following  wave-op- 
tlcail  argument. 

The  light  paaalng  thi'ough  a  grating  allt  of 
width  d/2  does  not  behave  like  a  ray  anymore 
alter  it  has  travelled  over  a  distance  Xj.  where 
the  diffraction  spread  XpX/(d/2)  equals  me  slit 
width  d/2  [eq.  (2)]. 

XpX/(d/2)  =  d/2  :  Xp  -  d2/4X  .  (2) 

Based  on  this  criterion  one  might  expect  Cp  to 
be  the  smallest  detectable  deflection  [eq.  (3)]. 

«P  =  d/2Xp  =  2X/a  .  (3) 

With  X  s  D.5  X  10*^  mm  and  d  s  10*^  this  leads 
to  which  is  only  moderately  good. 

Irortunxtely  the  method  is  capable  of  detecting 
angles  much  smaller  than  Cp  [eq.  (3)].  This  is 
due  to  the  Talbot  effect  [6],  which  is  also  known 
as  "Fourier  imaging"  [7]  and  "seif-lmaglng"  [8]. 
Talbot  discovered  about  135  years  ago  that 
images  can  be  formed  without  any  lenses  or 
mirrors  if  the  object  is  a  grating  which  is  il¬ 
luminated  in  collimated  monochromatic  light. 
These  "Talbot  images*  occur  at  distances  2d^A« 
id^/y,  etc.  behind  the  grating.  When  a  second 
grating  at  a  slight  angular  rotation  to  G1  is 
placed  into  the  plane  of  a  Talbot  image,  moird 
fringes  of  high  contrast  are  observed.  These 
moird  fringes  will  be  deformed  if  a  refractive 
phase  object  is  placed  for  example  close  to  the 
first  grating.  As  before  these  moird  fringes  in¬ 
dicate  lines  of  equal  deflection  by  the  object 
[eq.  (1)]. 

The  largest  possible  distance  x.j.  (or  maybe 
2X.J.)  depends  on  the  finite  width  B  *Nd  at  the 
first  grating.  At  Xf  the  first  grating  diffraction 
orders  have  moved  to  both  sides  by  half  of  the 
grating  width  B  [fig.  2;  eq.  (4)]. 

XjX/d  =  B/2  =  Nd/2  ;  Xx  =  M®/2X  .  (4) 

When  the  two  gratings  are  separated  by  xx  the 
smallest  detectable  deflection  Cx  is  now  smaller 
by  a  factor  ZN,  where  N  is  the  number  of  periods 
in  the  first  grating  [eq.  (5)]. 

Cx  =  <f/2xx  =  X/XW  =  X/B  =  Cp/2Af .  (5) 

We  have  performed  some  experiments  with 
this  "Talbot  interferometer"  as  shown  in  fig.  1. 

In  fig.  3a  the  object  is  a  flamo,  and  again  in 


fig.  3b,  but  now  with  the  second  grating  slightly 
rotated  around  the  optical  axis.  In  some  addi¬ 
tional  experiments  we  have  placed  a  telecentric 
lens  system  behind  the  Talb^  interferometer. 
When  introducing  a  spatial  filter  into  this  tele- 
centric  system  and  rejecting  everything  but  the 
zeroth  diffraction  order  from  the  two-grating 
system  one  observes  in  essence  the  second 
derivative  d^u/ax^  the  object  u(x,y).  When 
shifting  the  spatial  fiber  to  one  of  the  first 
grating  diffraction  orders  the  image  represents 
shearing  interferences  M(r -fXx/d,  y)  -u(;c,  y). 

In  a  white  light  one  obtains  beautiful  color 
fringes  which  are  unlike  ordinary  interference 
fringes.  For  example  for  a  specific  grating 
distance  the  image  contains  many  blue -orange 
fringes.  These  color  fringes  are  of  such  high 
contrast  that  E.  Lau  found  a  similar  setup  useful 
in  designing  tablecloths  [9].  Our  method  is  also 
applicable  for  the  detection  of  small  differences 
between  two  quite  irregular  but  similar  objects. 
One  begins  by  recording  photographically  the 
fringe  from  the  first  object.  After  development 
the  photograph  is  placed  where  the  fringes  had 
been  observ^.  The  first  object  is  now  replaced 
by  the  second  object,  which  might  actually  be 
the  first  object  but  somewhat  deformed.  The 
moird  fringes  between  subsequently  produced 
Talbot  interference  fringes  will  reveal  small 
differences  between  the  two  objects.  The  quanti¬ 
tative  evaluation  is  similar  to  that  for  Lau's 
diq>llgram  method  [10]  and  for  life-fringe  holo¬ 
graphic  interferometry. 
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There  are  two  coonon  aethoda  of  determining  the  degree  of  collimation: 

autocollimation  and  shearing  interferometry.  The  simplest  method  is  the  auto- 

collimation  technique,  but  it  only  Indicates  collimation  by  comparing  the  size 

of  the  source  with  its  image.  On  the  other  hand  the  spacing  of  the  fringes  in 

shearing  interferometry  is  a  direct  measure  of  the  degree  of  collimation.^ 

Recently,  Langenbeck  derived  a  method  that  used  two  tiny  corner-cube  reflectors 

to  sample  the  beam.  In  this  fashion  he  was  able  to  translate  the  measurement 

2 

from  fringe  spacing  to  fringe  rotation,  resulting  in  Increased  sensitivity. 

3 

The  method  presented  here  is  an  application  of  the  Talbot  Interferometer  and 
has  the  same  limitations  as  the  other  interferometric  techniques.  The  main 
advantages  of  this  technique  are  the  inexpensive  components  (two  Ronchl  rulings) 
and  the  relative  Insensitivity  to  component  alignment. 

As  Talbot  observed  in  1837,  an  image  of  a  grating  appears  at  Integral 
multiples  of  the  distance  2d^/X  (where  d  is  the  grating  period)  when  the  grating 
is  illuminated  by  a  plane  monochromatic  wave.  A  Ronchi  ruling  would  therefore 
be  Imaged  at  2dV^>  Ad^/X,  etc.  If  we  place  another  identical  Ronchl  ruling 
in  one  of  these  self-image  planes  (see  Figure  1)  moir8  fringes  will  be  formed, 

4 

as  observed  by  J.  Burch.  The  fringe  spacing  P  can  be  calculated  by  using  the 
vector  diagrams  described  by  Rogers^  to  be 

P  -  d/2  oin(0/2)  (1) 

where  6  is  the  angle  with  which  the  two  gratings  meet.  As  the  grating  62  is 
rotated  around  the  optical  axis,  the  fringe  spacing  increases  until  we  have 
uniform  brightness  (or  darkness)  when  grating  bars  of  61  and  62  are  parallel 
to  each  other. 

The  plane  wave  that  illuminates  61  of  the  interferometer  is  obtained  by 
placing  a  point  source  at  the  focus  of  a  lens  as  shown  in  Figure  1.  When  this 
lens  is  defocused  the  plane  monochromatic  wave  becomes  spherical.  Cowley  and 


1 


Moodle  have  shown  that  the  positions  of  the  self-image  planes  of  a  grating 
Illuminated  by  a  spherical  wave  are  given  by 

1/Zj  +  1/Zy  -  l/2JZx.  (2) 

vdiere  J  Is  a  positive  Integer,  1/Zy  Is  the  curvature  of  the  Incident  wave  at 
Gl,  Zj  Is  the  location  of  the  self-image  as  measured  from  G1  and  Zj  Is  the 
distance  2d^/A.  The  self-image  will  be  magnified  by  a  factor 

M  -  1  +  Zj/Zy  .  (3) 

This  Image,  when  superimposed  with  the  second  grating,  again  produces 
molrd  fringes.  Now  the  grating  G2  Is  rotated  around  the  optical  axis  until 
the  grating  bars  of  Gl  and  G2  are  parallel  to  laach  other.  In  this  case,  the 
Rogers  vector  diagram  gives  a  fringe  spacing 

P  -  did2/|di  -  dil 

where  d^  Is  the  period  of  the  self-image,  d2  the  period  of  the  second  grat.^ng, 
and  di  -  Md.  The  fringe  spacing  In  terms  of  the  grating  period  d  Is  therefore 

P  -  Md/|M  -  ij . 

Thus  the  frluge  spacing  Is  a  quantitative  measure  of  the  degree  of  colllmatlon 
since  |m  -  l|  Is  proportional  to  the  curvature  1/Zy. 

Now  we  want  to  discuss  the  accuracy  obtainable,  that  Is,  the  smallest 
detectable  deviation  from  perfect  colllmatlon.  If  the  colllmatlon  Is  perfect, 
the  magnification  M  of  the  Talbot  Image  Is  1,  and  hence  the  molr£  period  P  Is 
Infinite.  Thus  our  task  Is  to  find  how  small  must  the  period  P  become  In 
order  to  be  detectable.  Experience  has  shown  that  the  presence  of  a  molr^. 
effect  Is  clearly  visible  If  at  least  one  half  of  the  molr£  period  P  falls 


2 


within  the  Halted  field  of  observation  with  width  B.  Thus  one  fringe  will 
be  detected  whenever  the  number  of  lines  between  grating  G2  and  the  Image  of 
G1  differ  by  one-half.  With  divergent  Illumination  we  count  N  lines  of  G2 
in  the  observation  fie^  •  oi  width  B  and  N  -  1/2  lines  of  the  image  of  Gl.  Thus 
the  magnification  M  required  for  fringe  detection  is  bounded  by 


M 


N 


N  -  1/2 

and  in  convergent  Illumination, 

N 


«  1  +  1/2N, 


M  < 


«  1  -  1/2N. 


N  +  1/2 

By  combining  equations  (2)  and  (3)  we  have 

M  -  (1  -  2J  Zt/Zw)“^  «  1  +  2J  Zj/Zy. 


(A) 


(5) 


(6) 


If  M  falls  outside  the  bounds  of  equ.  (A)  and  (5),  no  fringes  will  appear, 
limiting  the  sensitivity  of  this  test.  Solving  equ. (6)  with  the  bounds  of 
equ.  (A)  and  (5)  we  find  that  the  minimum  detectable  field  curvature  1/Zy  Is 

ll/Zyl  ?.  (ANJ  Zx)“l 


where  the  absolute  signs  obviate  the  necessity  of  assigning  a  sign  convention 
for  field  curvature.  By  substituting  equ.  (2)  Into  the  above  equation  we 
obtain,  for  N  »  1, 

ll/Zyl  ^  <2N  Zj)-1  (7) 

From  equ.  (7)  It  follows  that  a  large  value  of  Zj  Is  desired  in  order  to 
detect  small  colllmatlon  defects  1/Zy  .  The  furthest  distance  to  which  we  can 
go  is  limited  by  the  "wallcoff"  of  the  first-order  grating  diffraction.  As 
indicated  in  Figure  2,  the  longest  distance  at  which  there  is  still  a  connected 
interference  field  of  width  D  is  at  Zj  ■  Dd/2X.  Perhaps  one  could  perform 
colllmatlon  tests  up  to  about  twice  this  distance,  but  we  will  restrict  ourselves 


3 


to  Zj  5  Dd/2X.  Inaertlng  this  Halt  Into  equ.  (7),  together  with  N  •  D/d,  we 
get 

jl/Zwl  i  2X/d2.  (8) 


The  curvature  of  the  wave  incident  on  G1  Is  related  to  the  defocuslng  of 
the  collimating  objective  by 


f  -  Zf 


If  -  ZqI 


(9) 
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Zg|f  -  Zol  +  Zof  Zof 

where  Zq  Is  the  distance  of  the  collimating  objectlvctfl||^oB  the  point  source, 

Zq  Is  the  distance  from  objective  to  Gl, and  f  Is  focal  length  of  the  objective 
lens.  The  approximation  Is  permissible  since  the  distance  Zq  <<Zof/|f  -  ZqI 
In  the  test.  Call  6  -  |f  -  ZqI  the  focusing  error,  and  substitute  equ.  (8) 
Into  equ.  (9).  We  have 

6  >  2A(f/D)2  (10) 


where  we  have  used  the  fact  Zq  ^  f.  Thus  the  focusing  error  Is  simply  related 
to  the  wavelength,  and  the  square  of  the  focal  length  to  beam  diameter  ratio. 

The  focusing  error  can  be  reduced  by  simply  rotating  one  grating  with 
respect  to  the  other,  thereby  producing  molr8  fringes  In  the  field  of  observa¬ 
tion.  In  this  fashion  we  change  our  detection  scheaie  from  the  condition  of 
uniform  brightness  to  a  system  cf  rotating  fringes.  It  Is  known  that  fringe 
detection  is  more  accurate  than  use  of  the  uniform  brightness  condition.  We 
rotate  both  gratings  In  opposite  directions  by  angles  6/2  with  the  y-axls.  If 
both  gratings  have  the  same  period,  molrd  fringes  will  appear  In  the  field  of 
obe^rvatlon  parallel  to  the  x-axis  with  spaclngs  given  by  equ.  (1).  When  the 
periods  are  unequal,  the  fringe  spacing  Is  given  by 

did2 

p  a  -  M  - - 

(df  +  d^  -  2did2  cose)^  |di  -  d2l 


A 


where  the  approzloation  is  for  sasll  angular  rotations  of  6,  and  dj,  d2  are  the 

grating  periods  of  G1  iaage  and  G2  respectively.  It  can  be  shown  by  using  the 

vector  arguaents  of  Rogers  that  these  fringes  of  unequal  grating  periods  are 

rotated  through  an  angle  ^  with  respect  to  the  z-azis,  and  the  direction  of 

this  rotation  depends  on  whether  the  llluaination  is  converging  or  diverging. 

Thus  the  ezperiaenter  knowtt  in  which  direction  to  aove  the  objective  lens.  This 

rotation  is  related  to  the  aagnification  M  by 

j-  _  cot6/2  tani 
cot6/2  -  tan^ 

Substituting  this  equation  into  equ.  (3)  we  obtain 


1  _  2  tand 

Zy  Zj  cot6/2  -  tan^ 


(11) 


Our  ability  to  aeasure  4  llwlta  the  sensitivity  of  this  test.  If  we  say  the 
minlaum  detectable  fringe  angle  4  Is  about  one-half  of  a  fringe  (see  Figure  3)* 
then  equation  (11)  reduces  to 


|l/Zy|  >  21/d2. 

This  lower  bound  is  identical  to  that  of  the  uniforw  brightness  condition, 
equ.  (8).  In  the  Isborstory  one-quarter  of  a  fringe  was  essily  discernible, 
producing  results  which  were  s  fsetor  of  two  better  than  predicted.  If  we 
want  to  achieve  even  higher  accuracy  we  way  use  photodetectors  instead  of  the 
eye.  Assualng  IZ  brlghtnecs  accuracy,  we  cen  detect  fringes  one  hundred  tlwes 
large?  than  the  field  of  view.  Hence  the  photoelectrical  approach  would 
improve  the  detectlbillty  of  colllaatlon  errors  by  a  factor  of  about  fifty 
over  that  which  was  predicted. 

The  visual  ezperiaents  have  confirmed  our  theoretical  performance  pre¬ 
dictions.  An  objective  lens  of  20  cm  focal  length  was  used  to  produce  a  5  cm 
beam.  The  light  from  a  helium-neon  laser  was  focused  with  a  3.9  mn  focal 
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length  lens  onto  a  <S  p  pinhole.  The  position  accuracy  was  0.001  cm,  whereas 
the  calculated  accuracy  was  0.002  cm.  This  Improvement  was  due  to  the  fact 
that  1/4  fringe  was  used  with  the  angular  measurement  technique. 

Last,  one  c^r.  -^mDrove  the  accuracy  of  the  test  by  increasing  the  distance 
Zj.  This  distance  was  limited  by  the  walkoff  as  shown  in  Figure  2.  By  the  addi¬ 
tion  of  one  or  two  mirrors  placed  idiere  the  borderline  rays  of  the  ^-'eroth 
diffraction  order  propagate  (Figure  2),  the  diffracted  first-order  beam  can 
be  reflected  back  into  the  zeroth-order  beam  to  create  the  shearing  inter¬ 
ferences.  This  should  Increase  Zj  and  the  accuracy  of  the  test  manyfold. 

However,  this  has  not  been  verified  experimentally. 

The  author  is  Indebted  to  Adolf  Lohmann  for  his  support  and  Interest  in 
this  work  and  for  his  helpful  suggestions  in  the  preparation  of  this  letter. 

This  work  was  performed  under  Air  Force  contract  AF-ri9628-69-C-0268 . 
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FIGURE  CAPTIONS 


1.  Talbot  Interferometer  setup  for  observing  fringes  due  to  defocusing  of 
collimator  objective  (Zq  ^  f ) .  Fringes  are  observed  on  G2. 

2.  The  self-imaging  effect  occurs  within  overlapping  orders.  Test  done  at 
distance  where  the  two  first  order  diffraction  spots  begin  to  separate. 

3.  Moird  fringes  which  are  rotated  through  angle  ^  because  of  defocusing  of 
collimator  objective. 
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